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We formulate a quantitative theory of non-local electron transport in three-terminal disordered 
ferromagnet-superconductor-ferromagnet structures. We demonstrate that magnetic effects have 
different implications: While strong exchange field suppresses disorder-induced electron interference 
in ferromagnetic electrodes, spin-sensitive electron scattering at superconductor-ferromagnet inter- 
faces can drive the total non-local conductance gi2 negative at sufficiently low energies. At higher 
energies magnetic effects become less important and the non-local resistance behaves similarly to 
the non-magnetic case. Our predictions can be directly tested in future experiments on non-local 
electron transport in hybrid FSF structures. 

PACS numbers: 74.45. -|-c, 72.25.Ba, 73.23.-b, 74.78.Na 



I. INTRODUCTION 

The phenomenon of Andreev reflection (AR)i is well 
known to be responsible for transport of subgap electrons 
across an interface between a normal metal (A'') and a su- 
perconductor (S). While this phenomenon is essentially 
local in hybrid proximity structures with only one A^S* in- 
terface, the situation in multiterminal devices with two or 
more NS interfaces (such as, e.g., NSN structures) can 
be more complicated because in addition to local AR elec- 
trons can suffer non-local or crossed Andreev reflection 
(CAR)^ . This phenomenon of CAR enables direct experi- 
mental demonstration of entanglement between electrons 
in spatially separated A^-electrodes and can strongly in- 
fluence non-local transport of electrons in hybrid NSN 
systems^.!* . 

Non-local electron transport in the presence of CAR 
was recently investigated both experimentallj^^— and 
theoreticalljj^i^iii"— demonstrating a rich variety of phys- 
ical processes involved in the problem. For instance, the 
effect of CAR on the subgap non-local conductance of 
NSN structures is exactly compensated by elastic cotun- 
ncling (EC) provided only the lowest order terms in NS 
interface transmissions are accounted for- . Taking into 
account higher order processes in barrier transmissions 
eliminates this feature and yields non-zero values of cross- 
conductanced.. One can also expect that interactions^^ or 
external ac biasi^ can lift the cancellation between EC 
and CAR contributions already in the lowest order in 
barrier transmissions. 

Another non-trivial issue is the effect of disorder. The- 
oretical analysis of CAR in different disordered NSN 
structures was carried out in Refs. [l3l - [l7l In particular, 
it was demonstratedii that an interplay between CAR, 
quantum interference of electrons and non-local charge 
imbalance dominates the behavior of diffusive NSN sys- 
tems being essential for quantitative interpretation of a 
number of experimental observations^^— . 



Yet another important property of both local and non- 
local Andreev reflection processes is that they essentially 
depend on spins of scattered electrons. Hence, CAR 
should be sensitive to magnetic properties of normal elec- 
trodes. This sensitivity was indeed demonstrated already 
in the first experiments on ferromagnet-superconductor- 
ferromagnet {FSF) structures^ where the dependence 
of non-local conductance on the polarization of ferro- 
magnetic terminals was found. Theoretical analysis of 
spin-resolved CAR was carried out in Ref. 3 i n the low- 
est order order in tunneling and in Refs. [l^ [l^ to all 
orders in the interface transmissions. This analysis re- 
vealed a number of non-trivial features of non-local spin- 
dependent electron transport which can be tested in fu- 
ture experiments. 

Note that previous work^ d^'^^ merely concentrated 
on ballistic electrodes whereas in realistic experiments 
one usually deals with diffusive hybrid FSF structures. 
Therefore it is highly desirable to formulate a theory 
which would adequately describe an interplay between 
disorder and spin-resolved CAR. This is the main goal 
of the present paper. The structure of our paper is as 
follows. In Sec. 2 we will formulate our model and 
outline our basic formalism of quasiclassical Green func- 
tions. This formalism will be employed in Sec. 3 where 
we present the solution of Usadel equations and derive 
general expressions for the non-local spin-dependent con- 
ductance and resistance for diffusive three-terminal FSF 
structures at different directions of interface magnetiza- 
tions. Concluding remarks are presented in Sec. 4 of our 
paper. 



II. MODEL AND BASIC FORMALISM 

Let us consider a three-terminal diffusive FSF struc- 
ture schematically shown in Fig. [TJ Two ferromagnetic 
terminals Fi and F2 with resistances r^-^^ and r^.^ and 
electric potentials Vi and V2 are connected to a supercon- 




FIG. 1: FSF structure under consideration. 



ducting electrode of length L with normal state (Drude) 
resistance tl and electric potential V^ = via tunnel bar- 
riers. The magnitude of the exchange field /11.2 = |^i,2| 
in both ferromagnets -Fi and F2 is assumed to be much 
bigger than the superconducting order parameter A of 
the 5'-terminal and, on the other hand, much smaller 
that the Fermi energy, i.e. A ^ /ii_2 <C cf- 

The latter condition allows to perform the analysis of 
our FSF system within the quasiclassical formalism of 
Usadel equations for the Green-Keldysh matrix functions 
G. In each of our metallic terminals these equations can 
be written in the fornt^S 



iDViGVG)^[n + eV.,G], G^ = 1, 



(1) 



where D is the diffusion constant, V is the electric poten- 
tial, G and Cl are 8x8 matrices in Keldysh-Nambu-spin 
space (denoted by check symbol) 
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e is the quasiparticle energy, A(r) is the superconducting 
order parameter which will be considered real in a super- 
conductor and zero in both ferromagnets, h = /ii(2) in 
the first (second) ferromagnetic terminal, h = outside 
these terminals and & — (fxi, (T2, (T3) are Pauli matrices 
in spin space. 

Retarded and advanced Green functions G^ and G^ 
have the following matrix structure 



G 



R,A 



qR.A 
_pR,A 



pR.A 
_qR,A 



(4) 



Here and below 2x2 matrices in spin space are denoted 
by hat symbol. 

Having obtained the expressions for the Green-Keldysh 
functions G one can easily evaluate the current density j 
in our system with the aid of the standard relation 



-— / Sp[T3(GVG)^]d£, 



(5) 



where a is the Drude conductivity of the corresponding 
metal and T3 is the Pauli matrix in Nambu space. 



In what follows it will be convenient for us to employ 
the so-called Larkin-Ovchinnikov parameterization of the 
Keldysh Green function 



G 



K 



G"f~fG^, f^h + rsh 



(6) 



where the distribution functions /^ and fr are 2x2 
matrices in the spin space. 

For the sake of simplicity we will assume that magne- 
tizations of both ferromagnets and the interfaces (see be- 
low) are coUinear. Within this approximation the Green 
functions and the matrix Cl are diagonal in the spin space 
and the diffusion-like equations for the distribution func- 
tion matrices /^ and /t take the form 

-DV (n^ir, £)V/T(r, e)) + 2E(r, e)fT{r, e) = 0, (7) 



where 



-i?v(i)^(r,e)V/L(r,e)) = 0, 



S(r,e) = -iAlmF'^, 
73^ = (ReG^y + (lmF«)', 



D 



^-^ReG^Y-fRei^^^' 



(8) 

(9) 
(10) 

(11) 



The function S(r, e) differs from zero only inside the su- 
perconductor. It accounts both for energy relaxation of 
quasiparticles and for their conversion to Cooper pairs 
due to Andreev reflection. The functions D^ and D^ ac- 
quire space and energy dependencies due to the presence 
of the superconducting wire and renormalize the diffusion 
coefficient D. 

The solution of Eqs. ([T])-® can be expressed in terms 
of the diffuson-like functions T)^ and 2?^ which obey the 
following equations 



-DV 



D^{r,e)W'D'^{r,r',e) 



-2±{r, e)'D^{r, r', e) =(5(r - r'). 



DV D^{r,e)V'D^ir,r',e) = 6{r - r') 



(12) 



(13) 



The solutions of Usadel equation ([T]) in each of the met- 
als should be matched at S'i^-interfaces by means of ap- 
propriate boundary conditions which account for electron 
tunneling between these terminals. The form of these 
boundary conditions essentially depends on the adopted 
model describing electron scattering at S'i^-interfaces. 
Here we stick to the model of the so-called spin-active 
interfaces^i which takes into account possibly different 
barrier transmissions for spin-up and spin-down elec- 
trons. This model was already extensively used for the- 
oretical description of different physical phenomena, in- 
cluding spin-resolved CAR in ballistic structures^^^ and 
Josephson effect with triplet pairin g^^i^^ . Here we em- 
ploy this model in the case of diffusive electrodes and 
also restrict our analysis to the case of tunnel barriers 



with channel transmissions much smaller than one. In 
this case the corresponding boundary conditions read^i 



+ ^[{^'mT3,G-},G+] +i^[&mT3,G+], 
+ -j^[{o-'mT3,G+},G-] +i-^[6-mr3,G_], 



(14) 



(15) 



where G_ and G_|_ are the Green-Keldysh functions from 
the left {x < 0) and from the right (a; > 0) side of the in- 
terface, A is the effective contact area, m is the unit vec- 
tor in the direction of the interface magnetization, a± are 
Drude conductivities of the left and right terminals and 
Gt is the spin-independent part of the interface conduc- 
tance. Along with Gt there also exists the spin-sensitive 
contribution to the interface conductance which is ac- 
counted for by the Gm-term, whereas the G^-teim arises 
due to different phase shifts acquired by scattered quasi- 
particles with opposite spin directions. 

Employing the above boundary conditions we can es- 
tablish the following linear relations between the distri- 
bution functions at both sides of the interface 



Aa+DldJ+T - Aa^D'^dJ-T 

= grif+T - f-r) + 9rnif+L - J^l), 

Aa+b^dJ+L = Aa-b'LdJ-L 

= hif+L - I-l) + 9mif+T - f-r), 



(16) 



(17) 



where gT, gL, and (jm are matrix interface conductances 
which depend on the retarded and advanced Green func- 
tions at the interface 



gr — Gt 
9l = Gt 



RcGI] (ReG 
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ReG2) (RcG^) - (RcF^) (RcFi^ 



G„6-m(ReG^) (rcG^) 



(18) 

r, 

(19) 

(20) 



The current density ^ can then be expressed in terms 
of the distribution function /t as 



4e 



Sp[b^VfT]de. 



(21) 



III. SPECTRAL CONDUCTANCES 

Let us now employ the above formalism in order to 
evaluate electric currents in our FSF device depicted in 
Fig. 1. The current across the first (SFi) interface can 



be written as 



/i = - / 311 (e) [/o(e + eVi) - /0(e)] de 



9i2ie)[Me + eV2)-Me)]de, (22) 



where /o(£) — tanh(e/2r), gn and 312 are local and non- 
local spectral electric conductances. Expression for the 
current across the second interface can be obtained from 
the above equation by interchanging the indices 1 <-> 2. 
Solving Eqs. dZ])-® with boundary conditions (fT6l) -(fT7 l) 
we express both local and nonlocal conductances 'gij{s) 
in terms of the interface conductances and the function 
T). The corresponding results read 



'■2m 



gnis) = {R^M'^+R^R^Rirn - R[RI 
512(e) = 321 (e) = {RT2M'' + RlRl2Rlm 

+ Rl2RlmR2m + ^12-^1 R2m)^, 

where we defined 



j(^T.L ^ flT^Lpr^L 



(nTX^2 
(-^12 j , 



t-^ = M^M'^+Rl^Rl^ - R^R^Rl 



'2Ri2Ri2RunR2v 



E>^ E>^ E>2 
Rl Rl R2r, 



(23) 
(24) 

(25) 
(26) 



and introduced the auxiliary resistance matrix 

RI = giT{£)[giTi£)giLie) - glmi^)]^^ 



<^\ 



OS 



(27) 



The resistance matrices i?J, R]( and R2 can be obtained 
by interchanging the indices 1 -H- 2 and T -(->■ L in Eq. 
(I27|) . The remaining resistance matrices i?]^2 and Rjm 
are defined as 



T,L, 



R 



12 



L _ f,T,L _ -PsPg' (ri,r2,£) 



^21 



(^s 



(28) 



Rjm = 5jm (e) [gjT {e)gjL (e) - g^n (e)] \ (29) 

where j = 1,2. The spectral conductance gtj can be 
recovered from the matrix gij simply by summing up over 
the spin states 



(e) = 2^P[5u(£) 



(30) 



It is worth pointing out that Eqs. ((23|) . (p4|) defin- 
ing respectively local and nonlocal spectral conductances 
are presented with excess accuracy. This is because the 
boundary conditions (|14[) -(|15 p employed here remain ap- 
plicable only in the tunneling limit and for weak spin 
dependent scattering |Gm|,|G<p| <C Gt- Hence, strictly 
speaking only the lowest order terms in Gmi 2 and G^p-^^ ^ 
need to be kept in our final results. 



In order to proceed it is necessary to evaluate the inter- 
face conductances as well as the matrix functions T>^ ^ s- 
Restricting ourselves to the second order in the interface 
transmissions we obtain 



girie) = Gr^vsivi^e) + G^^ — ^^ _ ^^ 

~ / ^ ^ " . ^ ^2 A20(£2 _ A2) 
giL [e] = Gti J^s (ri , e) - G^^ — ^TZT^ — 

5im(£) = Gmii>siri,s)&mi, 



Uiie), (31) 

Ui{s), (32) 
(33) 



and analogous expressions for the interface conductances 
of the second interface. The matrix function 

f7i(e) = ^|Re[Ci(ri,ri,2/i+)+Ci(ri,ri,2/ir)] 

Re [Ci(ri, n, 2/i+) - Ci(ri, n, 2/^^)] } (34) 



crmi 



with h^ = hi zL e defines the correction due to the prox- 
imity effect in the normal metal. 

Taking into account the first order corrections in the in- 
terface transmissions one can derive the density of states 
inside the superconductor in the following form 



i>sir,£) 



Ds A2 



as A2 



E 



GT.ReCs(r,n,2w«) 



i=l,2 



(35) 



where 



U! 



Ve2-A2, 

iV'A2-£2, 



£ > A, 
|£|<A, 



(36) 



Ve^ - A2, £ < A, 



and the Cooperon Cj{r,r',e) represents the solution of 
the equation 



{-DV^ - ie) C{r, r', £) = 5{r - r') 



(37) 



in the normal metal leads (j = 1,2) and the supercon- 
ductor {j = S). In the quasi-one-dimensional geometry 
the corresponding solutions take the form 



^ tanh (kjLj) 

Lj [Xj ,Xj,£) = Q r\ 1, ' J — 1 J ^: 

DjUjKj 



(38) 



smh fcs i - x ) smh fcsx , 

Cs{x,x,£)^ — . , . , X >x, (39) 

ksbsDsiiinh{ksL) 

where 55^1^2 are the wire cross sections and ^i 2,5 = 

yJ-ie/Di^2,s- 

Substituting Eq. ([35]) into Eqs. dSU and ([32]) and 
comparing the terms oc G'^ we observe that the tunnel- 
ing correction to the density of states dominates over the 



terms proportional to Ui which contain an extra small 
factor ^JKjh <c 1. Hence, the latter terms in Eqs. (|3T|) 
and ([32l) can be safely neglected. In addition, in Eq. (|35|) 
we also neglect small tunneling corrections to the super- 
conducting density of states at energies exceeding the 
superconducting gap A. Within this approximation the 
density of states inside the superconducting wire becomes 
spin-independent i>5(r,£) = a^vsir^e). It can then be 
written as 



vs(r,e) 



^|£2-A2 



+ 



i?S A20(A2 - £2) ^ 



1=1,2 



A2) 



GT,ReCs(r,r„2w^). (40) 



as A2 - £ 
Accordingly, the interface conductances take the form 

5iT(e) = g\L{e) = GT^vs(ri,e), (41) 

g\m^{e)=GmiVs{r-Y,e)crmx. (42) 

In the limit of strong exchange fields h\^i ^ A and 
small interface transmissions considered here the prox- 
imity effect in the ferromagnets remains weak and can 
be neglected. Hence, the functions T)^' (ri,ri,£) and 
I?2 ' (''2;'"2:£) can be approximated by their normal 
state values 

pP(ri,ri,£)=airA,J/i?i, (43) 

23^''^(r-2,r2,e)=a2r^,i/2?2, (44) 

rjv, =i,/(a,^,), j = l,2, (45) 

where r^q^ and rjVa are the normal state resistances of 
ferromagnetic terminals. In the the superconducting 
region an effective expansion parameter is Gria'^Csls)' 
where »'5s(^) = ^s(£)/(o'S'5's) is the Drude resistance 
of the superconducting wire segment of length ^5'(£) — 
^Dsl'l\i^^\. In the limit 

GT,,.r5,(£)«l, (46) 

which is typically well satisfied for realistic system pa- 
rameters, it suffices to evaluate the function 'D'g{x,x',e) 
for impenetrable interfaces. In this case we find 



A2 - £2 



T>Ux,x',e)^{^,^\, 



Cs{x,x',2u'^), \e\<A, 



Cs{x,x\0), 



(47) 



|£|>A. 



We note that special care should be taken while calculat- 
ing 'Dg{x,x',£) at subgap energies, since the coefficient 
D^ in Eq. (|S]) tends to zero deep inside the supercon- 
ductor. Accordingly, the function T>g{x,x' ,e) becomes 
singular in this case. Nevertheless, the combinations 
Rj{M^)~^ and Rf2{M^)~^ remain finite also in this 
limit. At subgap energies we obtain 



1 



■tni + ■rN2 + 



2(A2-£2)e' 



2\„d/(,s{e) ' 



(48) 



A2r5,(£)GT,GT. 



where d = |a;2 — a^il is the distance between two SF 
contacts. Substituting the above relations into Eq. ([M)) 

I 



we arrive at the final result for the non-local spectral 
conductance of our device at subgap energies 



512(e) = 521(e) 



A2-£2 



r^s{e)cxp[-d/£,s{e)] 



A2 2[rjvi + l/5Ti(e)][rjv, + l/5T2(e) 

Gml Gm2 A 



1 + mim2 



5Ti(e)5T2(e) A2 -£2 ^ e2 rN^+VN^ 



1 



A2 



r5<,(£)GTiGT2 exp[-d/^5(£)] 



|£| < A. (49) 



Eq. p9|) represents the central result of our paper. 
It consists of two different contributions. The first of 
them is independent of the interface polarizations mi_2. 
This term represents direct generalization of the resultii 
in two different aspects. Firstly, the analysiaii was car- 
ried out under the assumption r^i 25Ti,2(e) ^ 1 which 
is abandoned here. Secondly (and more importantly), 
sufficiently large exchange fields /ii_2 ^ A of ferromag- 
netic electrodes suppress disorder-induced electron inter- 
ference in these electrodes and, hence, eliminate the cor- 
responding zero-bias anomaly both in local^^^— and non- 
local^ spectral conductances. In this case with sufficient 
accuracy one can set gTi{£) = GTiVs(xi,e) implying that 
at subgap energies gxiis) is entirely determined by the 
second term in Eq. PO)) which yields in the case of quasi- 
one-dimensional electrodes 



2.5 



5Ti(e) 
5T2(e) 



A^GT,res(e) 

2(A2-£2) 

A^GT.rgs(e) 

2(A2-£2) 



Gt, +GT^e-'^^^'^''^ 



(50) 
(51) 



Note, that if the exchange field hi^2 in both normal 
electrodes is reduced well below A and eventually is 
set equal to zero, the term containing Ui{e) in Eqs. 
([5T|). ([5^ becomes important and should be taken into 
account. In this case we again recover the zero-bias 
anomalj«2^""— gnis) oc l/\/e and from the first term in 
Eq. (|49l) we reproduce the results^^ derived in the limit 
hi,2 -^ 0. 

The second term in Eq. (HH)) is proportional to the 
product mi'm2GmiG,n2 and describes non-local magne- 
toconductance effect in our system emerging due to spin- 
sensitive electron scattering at SF interfaces. It is impor- 
tant that - despite the strong inequality \G,ni\ ^ Gn ^ 
both terms in Eq. (J49| can be of the same order, i.e. the 
second (magnetic) contribution can significantly modify 
the non-local conductance of our device. 

In the limit of large interface resistances 
'^JVi 25Ti,2(e) <C 1 the formula P^ reduces to a 
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FIG. 2: (Color online) Local (long-dashed line) and non- 
local (short-dashed and solid lines) spectral conductances nor- 
malized to its normal state values. Here we choose rjvj = 
rN.2 = 5r5g(0), xi = L - X2 = 5^s{0), X2 - xi = Cs(0), 
Gt, = Gt2 = 4G,„i = 4G„2 = 0.2/r5g(0). Energy de- 
pendence of non-local conductance is displayed for parallel 
(P) miTn2 = 1 and antiparallel (AP) mim2 — —1 interface 
magnetizations. Inset: The same in the limit of low energies. 



much simpler one 

512(e) = 521(e) = ^ exp[-rf/^s(£)] 



r A2 



A2-£= 



5ri(e)5T2(e) + mim2GmiG„ 



A2 



(52) 

Interestingly, Eq. (I52p remains applicable for arbitrary 
values of the angle between interface polarizations mi 
and m2 and strongly resembles the analogous result for 
the non-local conductance in ballistic FSF systems (cf., 
e.g., Eq. (77) in Ref. [3). The first term in the square 
brackets in Eq. ([5^ describes the fourth order contribu- 
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FIG. 3: (Color online) Non-local resistance (normalized to 
its normal state value) versus temperature (normalized to the 
superconducting critical temperature Tc) for parallel (P) and 
antiparallel (AP) interface magnetizations. The parameters 
are the same as in Fig. [S] 
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FIG. 4: (Color online) The same as in Fig. 3 for the following 
parameter values: vmi ~ rN^ = 5rj3(0), x\ — L — X2 = 
5^5(0), a;2 - a;i = ^^(O), Gt, = Gt^ = 50G,„i = SOG™^ = 
0.025/r5^(0). 



tion in the interface transmissions which remains nonzero 
also in the limit of the nonferromagnetic leadaii. In con- 
trast, the second term is proportional to the product of 
transmissions of both interfaces, i.e. only to the second 
order in barrier transmission a^'^^ . This term vanishes 
identically provided at least one of the interfaces is spin- 
isotropic. 

Contrary to the non-local conductance at subgap en- 
ergies, both local conductance (at all energies) and non- 
local spectral conductance at energies above the super- 
conducting gap are only weakly affected by magnetic ef- 
fects. Neglecting small corrections due to Gm term in 
the boundary conditions we obtain 

gn{e)^RjiM^)-\ §22(6) ^ RJ{M^)-\ (53) 
gi2{e) ^ g2i{e) ^ RLiM^r', kl > A. (54) 

Eqs. ([53]) and (|54|l together with the above expressions 
for the non-local subgap conductance enable one to re- 
cover both local and non-local spectral conductances of 
our system at all energies. Typical energy dependencies 
for both 511(e) and (712(e) are displayed in Fig. [21 For in- 
stance, we observe that at subgap energies the non-local 
conductance gi2 changes its sign being positive for par- 
allel and negative for antiparallel interface polarizations. 

Having established the spectral conductance matrix 
gij (e) one can easily recover the complete I — V curves 
for our hybrid FSF structure. In the limit of low bias 
voltages these I — V characteristics become linear, i.e. 



h = Gn{T)Vi + Gi2{T)V2, 
I2 ^ G2i{T)Vi + G22{T)V2, 



(55) 
(56) 



where Gij{T) represent the linear conductance matrix 
defined as 



(57) 



^^ J cosh^ — 

2T 



It may also be convenient to invert the relations (|55p -([56 |) 
thus expressing induced voltages 14.2 in terms of injected 
currents h 2- 



Vi=RniT)h-Ri2iT)l2, 

V2 = -R2l{T)h+R22{T)l2, 

where the coefficients Rij{T) define local {i 
nonlocal [i 7^ j) resistances 



Rii[T) = 
Ri2{T) = R2i{T) 



G22{T) 



Gll(r)G22(T)-^12 

Gi2(T) 



GUTY 



Gii{T)G22{T) 



GUT) 



(58) 
(59) 

i) and 



(60) 
(61) 



and similarly for R22{T). In non- ferromagnetic NSN 
structures the low temperature non-local resistance 
Ri2{T — ?► 0) turns out to be independent of both the 
interface conductances and the parameters of the normal 
leads^^. However, this universality of R12 does not hold 
anymore provided non-magnetic normal metal leads are 
substituted by ferromagnets. Non-local linear resistance 
i?i2 of our FSF structure is displayed in Figs. El S] as 
a function of temperature for parallel (mim2 = 1) and 
antiparallel (mim2 ~ —1) interface magnetizations. In 
Fig. Owe show typical temperature behavior of the non- 
local resistance for sufficiently transparent interfaces. For 



both mutual interface magnetizations R12 first decreases 
with temperature below Tc similarly to the non-magnetic 
case. However, at lower T important differences occur: 
While in the case of parallel magnetizations R12 always 
remains positive and even shows a noticeable upturn at 
sufficiently low T, the non-local resistance for antiparal- 
lel magnetizations keeps monotonously decreasing with T 
and may become negative in the low temperature limit. 
In the limit of very low interface transmissions the tem- 
perature dependence of the non-local resistance exhibits 
a well pronounced charge imbalance peak (see Fig. H]) 
which physics is similar to that analyzed in the case of 
non-ferromagnetic NSN structures^ i^^'^'^ . Let us point 
out that the above behavior of the non-local resistance 
is qualitatively consistent with available experimental 
observations^. 



IV. CONCLUDING REMARKS 

In this paper we developed a quantitative theory 
of non-local electron transport in three-terminal hy- 
brid ferromagnet-superconductor-ferromagnet structures 
in the presence of disorder in the electrodes. Within 
our model transfer of electrons across SF interfaces is 
described in the tunneling limit and magnetic proper- 
ties of the system are accounted for by introducing (i) 
exchange fields ft.1^2 in both normal metal electrodes 
and (m) magnetizations mi 2 of both SF interfaces (the 
model of spin-active interfaces). The two ingredients 



(i) and (ii) of our model are in general independent 
from each other and have different physical implications. 
While the role of (comparatively large) exchange fields 
hi2 ^ A is merely to suppress disorder-induced inter- 
ference of electrons2^^— penetrating from a superconduc- 
tor into ferromagnetic electrodes, spin-sensitive electron 
scattering at SF interfaces yields an extra contribution 
to the non-local conductance which essentially depends 
on relative orientations of the interface magnetizations. 
For anti-parallel magnetizations the total non-local con- 
ductance gi2 and resistance R12 can turn negative at suf- 
ficiently low energies/temperatures. At higher tempera- 
tures the difference between the values of R12 evaluated 
for parallel and anti-parallel magnetizations becomes less 
important. At such temperatures the non-local resistance 
behaves similarly to the non-magnetic case demonstrat- 
ing, e.g., a well-pronounced charge imbalance peak23, in 
the limit of low interface transmissions. 

We believe that our predictions can be directly used for 
quantitative analysis of future experiments on non-local 
electron transport in hybrid FSF structures. 
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